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ABSTRACT 



Context. The weak field limit for a pointlike source of a f(R) oc /? 3/2 -gravity model is studied. 

Aims. We aim to show the viability of such a model as a valid alternative to GR + dark matter at Galactic and local scales. 
Methods. Without considering dark matter, within the weak field approximation, we find general exact solutions for gravity with 
standard matter, and apply them to some astrophysical scales, recovering the consistency of the same /(R)-gravity model with cos- 
mological results. 

Results. In particular, we show that it is possible to obtain flat rotation curves for galaxies, [and consistency with] Solar System tests, 
as in the so-called "Chameleon Approach". In fact, the peripheral velocity Vc is shown to be expressed as Voo = A yM, so that the 
Tully-Fisher relation is recovered. 

Conclusions. The results point out the possibility of achieving alternative theories of gravity in which exotic ingredients like dark mat- 
ter and dark energy are not necessary, while their coarse-grained astrophysical and cosmological effects can be related to a geometric 
origin. 

Key words. /(7?)-gravity - dark matter - gravity tests 
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1. Introduction 



Alternative th eories of gravity (Peebl es and Ratra (|2003); 
Padmanabhan d2003l) ; Copeland et al. (|2006)) are increasing 
as possible suitable alternatives to dark energy and dark mat- 
ter. Although the ACDM model is affected by many theoretical 
shortcomings (Carroll et al. dl992l) ). and, in general, dark en- 
ergy models are mainly based on the implicit assumption that 
Einstein's General Relativity (GR) is the correct theory of grav- 
ity. But the validity of GR both on large astrophysical and cos- 
molo gical s cales still remains to be accurately tested (see e.g. 
Will (2006)), and there is still enough room to propose dif- 
ferent theoretical schemes. A minimal alternative choice could 
be to take into account generic functions f(R) of the Ricci 
scalar curvature R. The task for these extended theories is to 
fit the astrophysical data without adding exotic dark ingredients 
(Klein ert and Schmidt (l2002h : Capoz ziello (I2002I) : Capozziello 
et al. d2003al) ; Odintso y and Nojiri d2003l) ; Cap ozziel lo et al. 
d2003bj); CaiTonetal. d2004 ; Allema ndi et al. d2004 ; Nojiri 
and Odintso y d2004 : Cognola et al. d2005l) : Capozziello and 
Francaviglia d2008l) ). 

In such a context, these higher order theories have obtained 
considerable attention in cosmology, since they seem to work 
well b oth in the late and i n the ea rly universe (see Capo zziell o 
d2002l) ; Capozziell o et al. d2003al) ; Odints oy and Nojiri d2003l) : 
Ca pozzie llo et al. d2003bl) : Carr oll et al. d2004 : Allem andi e t 
al. d2004 : Nojiri and Odints oy d2004 : Cognola et al. d2005h : 
Capozziello and Francaviglia d2008l) ). It is also possible to show 
that f(R) theories can play a major role at astrophysical scales, 
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due to the modifications of the gravitational potential in the 
low energy limit. Such a corrected potential reduces to the 
Newtonian one at the Solar System scale and could also offer 
the possibility of fitting galaxy rotation curves and galaxy clus- 
ter potentials without the need for large amounts of dark matte r 
(Capozziello et al. (2004); Milgrom dl983l) . Be kenste in d2004 : 
Capoz ziello et al. d2006l) ; Capozziello et al d2007l) ; Sobouti 
d2007l) ; Frigeri o Mar tins and Salucci d2007l); M endoza and 
Rosas-Guevara d2007l) ; Capozziello et al. (2008a)). However, 
extending the gravitational Lagrangian may give rise t o man y 
problems. These theories may have instabilities (Faraoni (2005); 
Cognola and Zer bini d2006l) ; Cognola et al. d2007l) ). ghost-like 
behavior (Stelle dl978l) ). and they still need to be matched with 
data from the low energy limit experiments that are well under- 
stood by GR. 

In summary, adopting /(/?)-gravity leads to interesting re- 
sults, first of all at cosmological and galactic scales, even if, up 
to now, it has not been possible to select only one theory (or class 
of theories) good at all scales. There has been much work on this 
(Capo zziello et al. d2005l) ; H u and Sawicki d2007l) ; Starobinsky 
(2007); Nojiri and Odintsov (2007)), but all the approaches are 
indeed phenomenological and are not based on a fundamental 
conservation or invariance principle of the theory. 

In this paper, we propose a specific expression of the func- 
tion f(R) of the Ricci scalar curvature R, nameljQ f(R) = — | — 
R\ 3/2 , which comes from the need for the existence of a Noether 
symmetry for f(R) cosm ological models (de Ritis et al. £[990); 
Capozziello et al d2008bl) ). The cosmological solutions of the 
Einstein field equations related to su ch a ch oice for f(R) have 
been analyzed in Capozziello et al. (2008bJ, and it turned out 



1 where the combination of minus signs is only due to our conven- 
tions, since we start from the metric (- + ++) and we obtain R < and 

G ef f = -l/f. 
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that this model admits a dust-dominated decelerated phase, be- 
fore a late time accelerated phase, as needed by the observational 
data. Here, we study the low energy limit of such a solution, in 
the case of a point-like source. We consider the Schwarzschild- 
like spherically symmetric metric in such a way that, in the weak 
field limit, the Newtonian potential is modified by adding a log- 
arithmic term. A similar treatment has been proposed in Sobouti 
(2007), where instead the starting point consists of introducing a 
specific hypothesis on the metric and thereby deducing the form 
of f(R) (resulting in a power-law), so fitting observational data 
on speeds of peripheral sta rs in s piral galaxies, as first reported 
in San ders and McGough (120021) and then selected in Sobouti 
(2007). Unfortunately, this procedure leads to a very peculiar 
choice of f(R). It contains parameters which must be adjusted to 
the mass of the gravitational source. Therefore, f(R) cannot be 
universal. 

He re, we find results that are very similar to those in Sobouti 
(2007) from the observational point of view, but do not exhibit 
the above problems. In this preliminary work, we only consider 
the weak field generated by a point-like source. Clearly, shrink- 
ing a whole galaxy (or cluster) to a point is a very crude ap- 
proximation. Our aim, therefore, is to show that the model can 
nonetheless work, without trying to obtain a strict correspon- 
dence with observations. 

The treatment of point-like source in f(R) theories is non- 
trivial and, in a sense, can be considered an ill posed problem. 
The reason is that we cannot disregard the properties of the ex- 
tended object that is generating the field. Unlike what happens in 
GR, the choice of the integration constants is not necessarily in- 
dependent of its peculiarities, like density, equation of state etc. 
We should therefore solve the equations for the inner metric and 
then match with the exterior. This procedure is valid for a star 
but much less meaningful for a galaxy. In any case, we do not 
expect to have a full prediction of their functional dependence. 
In fact, also in standard GR the linear dependence on mass of 
the constant is obtained only from the observational link with 
Newtonian gravity. This link with observations is precisely what 
is lacking in the case of stars, and is only preliminarily stud- 
ied here. Therefore, in the following we retain the assumption 
of point-like source, dedicating Sec. 6 to a deeper discussion on 
this point. 

There is another limit of the analysis presented here, lying in 
the weak field approximation assumed from the beginning. The 
asymptotic Minkowskyian behavior cannot recovered and, most 
of all, does not shed light on the singularities of the metric. It 
is therefore not possible to say what are the modifications of a 
black hole so generated. 

In Sect. 2, we work out the basic equations of our model, and 
in Sect. 3 we study peripheral speeds in spirals. In Sect. 4, we 
discuss some tests in the Solar System and, in Sect. 5, we com- 
ment on gravitational lensing and MOND. In Sect. 6 the connec- 
tion between the law for a pontlike source and an extended body 
is briefly discussed. In Sect. 7, we give our conclusions. 



2. Weak field approximation 

Our theory of gravity is determined by the Action: 



due to the different starting equations. We thus set: 



(1) 



where R is the Ricci scalar and we define f(R) = — | - R^ 2 . 

As said above, we consider a static, spherically symmetric 
metric, which will differ from the standard Schwarzschild form, 



ds l 



B(r)dt 2 - A(r)dr 2 - r 2 (d6 2 - sin 2 0dcf> 2 ) , (2) 



with A(r) and B(r) radial functions to be determined through the 
modified Einstein field equations. 

It is important to observe that, as the angular coordinates are 
dimensionless, we also use a dimensionless distance, so that we 
have: 

r = p/r s , (3) 

where p is the physical radius and r s is a suitable scale. The 
choice of this scale is a delicate point: we could decide to use a 
universal scale or one that is specific for the particular situation. 
This second choice does not affect the universal character of our 
Action, but the situation will become clearer in the following. 

Proceeding in this way, all quantities are dimensionless, in- 
cluding R and the result for velocities. In order to restore the 
appropriate dimensions, we should therefore multiply by an ap- 
propriate fundamental Action, i.e. some numerical multiple of 
h. As we are in vacuum, this clearly does not affect the equa- 
tions and their solution. On the other hand, this constant is of 
course not irrelevant and has an influence on the coupling with 
the test particle (peripheral star or other object). In our case this 
arbitrariness will be resolved by restoring the physical units in 
an appropriate way when defining the observational objects and 
fixing the constant C3 (see below), with the prescription that we 
should obtain ordinary Newtonian gravity at small scales. 

In a weak field, of course, the functions A(r) and B(r) are 
practically identified by means of their corrections to what we 
expect in the Newtonian case for astrophysical situations. First 
trying to understand how we can modify the Newtonian poten- 
tial, we write the first function as B(r) = 1 — 2ey(r)/r + 0(e 2 ), 
where e is a suitable small parameter. Analogously, we also as- 
sume A(r) — 1 + 2ex(r)/r + 0(e 2 ). Being in a weak field, r 
cannot extend to infinity, but we must have 2ey(r)/r <k 1 and 
2ex(r)/r <sc 1. 

We can obtain our results depending on r, as for instance: 



R(r) = 



e(4x'(r) + 2ry"(r)) 



+ 0(e 2 ), 



(4) 



where prime denotes a derivative with respect to r. It is also: 
Rtt{r) = ^l + 0^). 



(5) 



To write equations, we have to vary the Action S with respect 
to the metric tensor, always remembering that we are studying an 
approximate situation. We also need to note that 



The master equation is: 

dR\gtt ) r 

that is: 

J-^i-M + Sry") 

°=~ F > 

V2r 2 

while the trace equation is given by: 

3aR'+R— -If = 0, 
dR J 



(6) 



(7) 



(8) 



(9) 
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which is equivalent to: 



(9 V^(4jc" 2 + y" 2 - 8x'x"' - 8x'y"' + r 2 y"' 2 + 4jc"(y" + ry'") + 
- Arx'y w - 2ry"(2x"' + y'" + ry (4) ))) X 



that the quantity under the square root is dimensionless, we ad- 
just the dimensions by means of the light speed c. This will fix 
the numerical values of the constants: 



4V2r 



4 / 2x' + ry' 



+ 0(e 3/2 ) = 



(10) 



Neglecting higher order terms and dividing the master equa- 
tion by e 3 '' 2 and the trace equation by e we finally find 



2x' + ry" 

^-(-8x' +5ry") =0, 

r L 

Ax" 2 + y" 2 - 8x'x"' - %x'y"' + r 2 y"' 2 + 4x"(y" + ry'")- 



(11) 



4rx'y m - 2ry"(2x" + y'" + ry {4> ) = . 



(4h 



(12) 



Such equations appear only difficult to handle, since Eq. (fTTb can 
be algebraically solved for x'(r) in terms of y"(r). Discarding the 
solution that removes the denominator in Eq. (2.10), we are then 
left with 

x'(r) = (5/8)ry"(r). (13) 

Substitution of this expression and its derivatives in Eq. ( TT2l 
leads to the general exact solution: 



1 



1 



y{r) = ^C 2 r L + C 3 + C 4 r- -C 2 r cos(2d)+ 
1 3 1 

— C 2 r cos(2Ci) + -C 2 rcos(2d)logr + 

1 1 3 

-iC 2 rsin(2Ci) + +^ lC 2 r sin(2Ci)- 

1 

-iC 2 rsin(2Ci)logr, 



(14) 



where C\, C 2 , C3, and C4 are arbitrary (complex) constants. If 
we limit ourselves to considering C 2 , C3, and C4 as real con- 
stants, we can understand that posing sin(2Ci) - ia and thus 
cos(2Ci) = Vl + a 2 (with a e R) makes y(r) real: 

y(r) = ^C 2 r 2 + ^C 2 (-a + Vl+a 2 )r 3 + C 3 + 



\r(-(C 2 a + C 2 Vl + a 2 - 4C 4 )+ 



C 2 (a + Vl + a 2 ) log r) . 



(15) 



Here, the constant C4 is multiplied by r and, since the potential 
is obtained from y(r)/r, it gives rise to a constant term. So, we 
can set C4 = in the following. On the other hand, the presence 
of a term giving the Newtonian potential explicitly depends on 
C3, which then needs to be nonzero. 

From now on, e will be incorporated into the integration con- 
stants without any loss of generality. 



3. Peripheral speeds 

We are now in a position to study how the speed of a test star 
behaves at the periphery of a pointlike galaxy in our model (i.e. 
at large distances from a point source of great mass), once it is 
subjected to the dimensionless potential <l>(r) = ^r. We write 
the speed v of a generic test star in this potential. Remembering 



dO 



c -C 2 r(a + Vl + a 2 ) + r(3 + r(-a + V 1 + a 2 )) + 12C 3 

d 



3r 



(16) 



The comparison with the Newtonian case is obtained when a = 
C 2 = 0. The value of C3 has to be C3 = ™. On the other hand, 
to get increasing velocities we need C 2 to be negative, and we 
define C 2 = ~vo 2 /3/c 2 . Defining Vo = y/l2GM/ r s , we can then 
rewrite v = v(r): 



v(r) = v 



' 1 + fir(3{a + Vl + a 2 ) + r{3 + r(-a + Vl + a 2 ))) 



(17) 

with a and /3 being dimensionless constants. Another (dimen- 
sional) constant is hidden in r — p/r s . As said in Sec. 1, the 
value of these constants should be linked to the inner peculiar- 
ities of the body. Here, we try a much simpler procedure, and, 
by means of suitable guesses, try to see if there is a convenient 
choice that reproduces the observations. The formulation of the 
solution so as to have dimensionless quantities was very helpful. 

The first assumption is a » 1, which clearly is arbitrary and 
can be justified only a posteriori. 

When a » 1, the expression under the square root reduces 
to (l/12)r~' + (l/4)y8r + (l/2)a/3. In fact, we find a correction 
of what is usually expected in the Newtonian case, since the first 
term is the Keplerian one, the third gives a constant speed, and 
the second must stay small. The reason is that, as said above, we 
must keep 2y(r)/r <K 1. This means that, even if a » 1, f3 must 
be small. However, it cannot simply be set to zero, for this would 
also remove the other correction term. 

The second assumption is therefore that /3, being so small 
with respect to a, we can safely neglect the increasing term. A 
rough estimate of the relevant correction term for a galaxy is: 



= vo a / —ctj3 — 23 1 ya^S km s 



(18) 



(where we have considered M = 10 10 M o , r s = 5Kpc, and G 
( j v ). indicating that y[afi a 1. 

Being dimensionless quantities, a and f3 can be universal 
numbers or arbitrary functions of GM/(c 2 r s ), i.e. the only num- 
ber we can form with the physical quantities at our disposal. But 
this is true only in the point-like source model. Another dimen- 
sionless quantity that can be considered for an extended body is 
of course its radius, again divided by r s , and many others. 

The third third assumption is thus that we consider a and 
/3 slowly varying with respect to these kinds of parameters. In 
other words, we expect a difference to appear only when a very 
large change of scale is considered, e.g. passing from the Galaxy 
to the Solar System, which is discussed in Sec. 6. 

The scale radius r s , having only the function of fixing the 
units, can be chosen arbitrarily and we set it as a multiple of 
the Schwarzschild radius, r s = KGM/c 2 , K being an unspecified 
positive (large) number (if K is kept fixed, then r s is chosen ac- 
cording to the first issue (see above); if, on the contrary, we want 
r s to be specific, we may suitably adjust the value of K. We shall 
see that this makes no difference). 
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We thus get: 



v 2 (p) 



6c 2 aB GM 3c 4 B 

1 1 : O , 

K p K 2 GM 



(19) 



that is, a constant term and a Keplerian one, plus one increas- 
ing contribution which can be neglected for sufficiently small /3 
and p. This shows that the relevant parameter is the combina- 
tion afi/K, while 3/3/K 2 must be negligible. We have thus re- 
duced our problem to the determination of a unique parameter 
C = a/3/K. We want however to stress that this simplification 
is due to the above assumptions. Moreover, it is clear that as- 
sumption three is more appropriately applied to C instead of its 
separate elements. 

Any attempt to modify the Newtonian law, on the Galactic 
scale, has to cope with the justification of the Tully-Fisher em- 
pirical relation. This obliges one to set the parameters of the 
modified part as depending on the mass of the source, which 
is impossible if one starts from a universal modification. On the 
contrary, it is appropriate here, where the parameters are spe- 
cific to the problem. This is indeed a very good feature of our 
approach. Moreover, we have complete freedom in the choice 
of the dependance. A reasonable assumption is to set C = 
proportional to some power of the mass M, made dimension- 
less by dividing by a suitable reference mass. Thus, we define 
AM" = 3af3/ K, where n is a pure number and M\ is referred to 
1O 1O M . (Such a definit ion of A is the same as the one for the 
parameter a in Sobouti d2007l) . 

If we neglect the term linear in p, we find: 



v 2 (p) 



c 2 M\A + 



GM 



(20) 



which coin cides w ith the relevant part of the expression found 
in Sobouti d2007l) . On the other hand, the constant term can be 
written as: 



= 21 1985 -M^ 2 VI fern j 



(21) 



A determination of A and n can be found by fitting this formula 
to the data on flat r otation curves from Sanders and McGough 
( 2002) as in Sobouti d2007l) . where a data list of 3 1 spirals is suit- 
ably selected and reported, including radii, total masses, asymp- 
totic orbital speeds, and velocity curves for each of those galax- 
ies. The fit is made with a non-linear regression algorithm and, 
as expected, we obtain the same result, as illustrated in Fig. 1. 
We obtain: 



« = 0.49 + 0.02 A - (3.1 + 0.1) ■ 10~ 
compatible with n = 0.5 at lcr. This yields: 

v oc M 1/4 , 



(22) 



(23) 



which is the phenomenological Tully-Fisher relation. We have 
thus proved that even in this crude model the flatness of the rota- 
tion curves can be obtained (a more detailed verification would 
require to treat the galaxy as an extended object); the immediate 
question to answer is now the validity of the theory at the Solar 
System level. 

4. Solar System tests 

4. 1 . Kepler's law 

We now use the result of the above section for in the Solar 
System, with the slight simplification n = 1/2. This is a ma- 
jor extrapolation, as nothing guarantees that a semi-empirical 



250 



200 



150 



100 



50 



o L 



10 



15 



20 



M 



Fig. 1. Fit of Eq. (Bil l versus data in Sanders and McGough 
(2002) as selected in Sobouti d2007l) . On the x-axis, masses of 
luminous matter in units of 10 10 M o , peripheral speeds in km s~ x 
in the y-axis. Error bars are not indicated as they are not reported 
in the reference papers. It is not crucial, since we are discussing 
only a preliminary estimate of the parameters. 



law like v (p) ~ (l/2)c M"A + GM/p can be used properly at 
such different scales, with the same values for the parameters. 
Possible arguments against this are: 

- The power law guess may be not the right one. It could in- 
deed be more complicated and show effects at small scales. 

- We do not know whether a » 1 at the Solar System level or 
not. 

- The parameters used for the Solar System tests are those used 
for galaxies, and they are already roughly estimated there. 

- The values used for G and M are the standard ones, which 
is not certain, since they are estimated within the Newtonian 
framework. 

- The orbits of the planets around the Sun are only approxi- 
mately circular. 

- The Solar System is a many body system, where perturba- 
tions should be accurately taken into account. (The observed 
deviations from Kepler's laws are in fact explained by the 
existence of perturbations induced by other planets.) 

- The passage from a pointlike source to an extended one is not 
trivial, as said above. This point seems the most important 
and is treated specifically below. 

If we use Eq. d20l i for the calculation of the Earth's orbital 
speed, we find VEanh - 29.8 kms~ x , with an error of ~ 1.4%, 
which is clearly large at this level. 

In the Solar System we can also test Kepler's third law 
T 2 /p 3 = Co (with T the period and Co a constant), where: 



2J1 



4n z p- 

Vplanet 2 



(24) 



and if we use v p i met = GM Q /p, then Co = 4n /(GM G ). When 
we instead use Eq. ( f20b for v p i anet 2 , we find a non-constant result: 



C(p) = 



An 1 



GM + 1.5 • 10- 7 c 2 Mi /2 p 



(25) 



The observed values of Co and those obtained from Eq.d25ll for 
Solar System planets are plotted in Fig. 2 against the semimajor 
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Sec. 5 this problem will be addressed by the introduction of a 
"Chameleon" mechanism. 
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Fig. 2. Comparison between the values of the observed Cos for 
the planets of the Solar System and what is obtained through Eq. 
d25l ). x-axis: the semi-major axes p of the planet orbits in the 
Solar System (in AU units); y-axis: the observed values of Co 
(points) and the values of C(p) computed from Eq. ( fZ5b (line). 

axes of the planet orbits. The deviation from observation is made 
evident by the absence of the downwards trend, which is instead 
obtained from Eq. d25l >. However, this is not a conclusive argu- 
ment, as the trend could be altered by the complicated effects of 
the many body system. 

4.2. PPN parameter 

In the study of modified (with respect to GR) approaches to grav- 
itational physics, it is usual to ask which are the new values for 
the PPN parameters, in particular of the most important and stud- 
ied y, which must be very near to 1 in the Solar System (and 
exactly 1 in GR). 

Since we have developed our weak field approximation only 
up to the first order, there is some doubt about the applicability 
of the formalism. We have a modification of the A(r) function in 
the metric via the function x(r) + 1, so that we may set: 



. . , 5 Ac 2 VMT 
y = x(r) = 1 — — p , 



4-GM 



finding: 



J earth 



0.999 . 



(26) 



(27) 



Thus, y turns out to be 1 up to 10~ 3 , which is indeed still far from 
the 10~ 5 approximation of the Cassini experiment. Moreover, we 
are not able to obtain other parameters, due to the linear devel- 
opment in e of the metric. 

The reason why we do not use the conformal transformation 
to pass to Einstein frame is that the equivalence can be lost in the 
case of the weak field approximation due to the fact that a f(R)- 
fourth-order theory has different gauge conditions in this limit 
with respect to a second-order Einstein theory plus a scalar field, 
which could lead to unphysical results. For a detailed discussion 
of this point see Capozziello et al d2007bl) . In order to avoid this, 
we adopted the Jordan frame for all calculations. 

We obtain a value for y which is not the same throughout 
the Solar System. On the other hand, the numerical codes that 
give an estimate of the PPN parameters are based on constant 
values. A precise computation should instead take this feature 
into account. We conclude that the transport of the correction 
as it is at the Solar System level seems to be unsatisfactory. In 



5. Other tests 

Still fixing, for simplicity, n = 1 /2, we also set K = 1 as it intro- 
duces a constant term in the potential, which is clearly irrelevant. 

5.1. Gravitational lensing 

Because it is intimately related to the underlying theory of grav- 
ity in its Einstein formulation, modifying the Lagrangian of the 
gravitational field also affects the theory of gravitational lensing. 
We therefore investigate how gravitational lensing works in the 
framework of higher order theories of gravity. On the one hand, 
one has to verify that the phenomenology of gravitational lens- 
ing is preserved in order to not contradict those observational 
results that do agree with the predictions of the standard theory 
of lensing. On the other hand, it is worth exploring whether devi- 
ations from classical results of the main lensing quantities could 
be detected and work as clear signatures of a modified theory 
of gravity. As a first step towards such an ambitious task, here 
we investigate how modifying gravity affects the gravitational 
lensing in the case of a point-like lens. 

Indeed, the basic assumption in deriving the lens equation is 
that the gravitational field is weak and stationary. In this case, 
the spacetime metric reads: 



</.v- = |l + ^W 



(28) 



where O is the gravitational potential and, as usual, we have ne- 
glected the gravitomagnetic term. Since light rays move along 
the geodetics of the metric in Eq. ( |28T i. the lens equation may 
be simply derived by solving ds 2 = 0. Such a derivation holds 
whatever the theory of gravity is, provided that one can still write 
Eq. ( l28l i in the approximation of weak and stationary fields. As 
a fundamental consequence, the lensing deflection angle will be 
given by the sam e form al expression found in general relativity 
(Schneider et al. dl992l) ; Petters et al. d200ll) ) 



a 



c 2 f 



v ± a>di, 



(29) 



with: 

V x = V - e(e-V) , (30) 

where e is the spatial vectoiQ tangent to the direction of the light 
ray and dl - yjSjjdx'dx' is the Euclidean line element. The in- 
tegral in Eq.d29l should be performed along the light ray tra- 
jectory which is a priori unknown. However, for weak gravita- 
tional fields and small deflection angles, one may integrate along 
the unperturbed direction. In this case, we may set the position 
along the light ray as: 

r = £ + le, (31) 

with £ orthogonal to the light ray. 

Assuming the potential only depending on r = \r\ = (£ 2 + 
l 2 ) 1 ! 2 (as for a point-like or a spherically symmetric lens), the 
deflection angle reduces to 



a= 2 4 r(— W 

c 2 J-oo \r dr j 



(32) 



2 Here, quantities in boldface are vectors, while the versor will be 
denoted by an over - hat. 
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where we have assumed that the geometric optics approxima- 
tion holds, the light rays are paraxial and propagate from infinite 
distance. Eq.d32b allows us to evaluate the deflection angle, pro- 
vided that the source mass distribution and the theory of gravity 
have been assigned, so that one may determine the gravitational 
potential. Here, we consider only the case of the point-like lens. 
Note that, although being the simplest one, the point-like lens is 
the standard approximation for stella r lense s in microlensing ap- 
plications (e.g. see Schneider et al. dl992l) ). Moreover, since in 
the weak field limit a is an additive quantity, the deflection angle 
for an extended lens may be computed integrating the point-like 
result weighted by the deflector surface mass distribution under 
the approximation of a thin lens (i.e., the mass distribution ex- 
tends over a scale that is far smaller than the distan ces between 
observer, lens and source) (Schneider et al. dl992l) ). Given the 
symmetry of the problem, it is clear that we may deal with the 
magnitude of the deflection angle and of the other quantities of 
interest rather than with vectors. 

In the approximation of small deflection angles, simple geo- 
metrical considerations allow us to write the lens equation as 



2ll 



(33) 



which gives the position 9 in the lens plane of the images of the 
source situated at the position 9 S on the source planfl Note that 
the lens and the source planes are defined as the planes orthog- 
onal to the optical axis, which is the line joining the observer 
and the center of the lens. Here and in the following, D/, D s , Di s 
are the observer - lens, observer - source, and lens - source angu- 
lar diameter distances, respectively. In order to evaluate the de- 
flection angle, we need an explicit expression for the gravita- 
tional potential <£ generated by a pointlike mass. In our case the 
modified dimensional Newtonian potential, dropping unimpor- 
tant constant terms and with the above assumptions, is: 



GM 
P 



<D = — - -2c 2 A yfmiog 



\GM 



(34) 



where we should remember that M\ is measured in units of 
10 10 M o . 



The first term yields, of course, the usual deviation: 

4GM 



S = 



£c 2 



(35) 



where £ is the angular impact radius. Computing the integral 

X+CO 
V ± <&dl leads to the additional deviation 
-co - 1 - 



Si = 2nA y/Mi , (36) 

which is independent of Now, in the case of the Sun we get 

6q ~ 8.5 x l(T 6 rad ; 6 X ~ 1.7 x 10 _1I rad, (37) 

so that the correction is irrelevant. On the other hand, in the case 
of a galaxy, assuming M\ — 1 and £ = lOKpc, we get 

S * 1.9 X l(T 6 rad ; Si * 1.7 X l(T 6 rad , (38) 

so that the correction now dominates and induces an overestima- 
tion of the mass, if the deviation is instead attributed in the usual 
wajQ. 



3 Both 6 and 6 S are measured in angular units and could be redefined 
as 9 = f/A and 6 S = r]/D s with f and 77 in linear units. 

4 Here, we limit ourselves to observing that, of course, the actual 
computation should be much more complicated, not only because of 
the fact that a galaxy is an extended object, but also because of the 
necessity to recompute the cosmological angular diameter distances 



The point-like lens equation differs from the standard gen- 
eral relativistic one for the second term d36l l. Should this term 
be negligible, all the usual results of gravitational lensing are re- 
covered. It is therefore interesting to investigate in detail how the 
corrective term affects the estimate of observable related quanti- 
ties since, should they be detectable, they could represent a sig- 
nature of R 1 ^ 2 gravity. Since we are considering the pointlike 
lens, a typical lensing system is represented by a compact ob- 
ject (both visible or not) in the Galaxy halo acting as a lens for 
the light rays coming from a stellar source in an external galaxy, 
like the Magellanic Clouds (LMC and SMC) or Andromeda. It is 
easy to show that, in such a configuration, the standard Einstein 
angle 9e.gr and the image angular separation are of the order 
of a few xlCT 5 arcsec, so that we are i n the r egime known as 
microlensing (see Mollerach and Roulet (120021) ). 

In the standard case the lens equation may be solved analyt- 
ically and one gets two images with positions given by 



(39) 



In the current case, the the lens equation may be conveniently 
written as 



D. 9, 



Mil?, 



(40) 



s Oe,GR 



with & - 9/9e,gr an d we nave defined 9e,gr as the Einstein an- 
gle, which in the general relativity case is given by 



\AGmDu 
' c 2 D,D s 



(41) 



We still get two images on the opposite sides of the lens, with 
one image lying inside and the other one outside the Einstein 
ring. The geometric configuration is therefore the same as in the 
standard case, but the positions are slightly changed. The images 
positions are given by 



i (A + f ± ^4 + (A + tf] 



(42) 



with A = In A 



Oh 



Be.gkD,, 



To quantify this effect, we evaluate the percentage deviation 
relative to the Einstein case. 9r = 1/2 (A + V4 + A 2 ), which dif- 
fers from 9e,gr- 

5.2. Equivalence with MONO 

Following Sobouti (120071) . we want now to show that our model 
is effectively equivalent to th e well known MOND theory (see 
Milgrom dl983l) : Bekenstein (2004)), which postulates a modi- 
fication of Newtonian dynamics: 



(43) 



F — mafi — , 



where p(x) is a suitable function, subject to the conditions: 

(i(x) = l(x»l) ; fx(x) -> (x -> 0) . (44) 

F is the usual gravitational force F = GM/p, while ao is a uni- 
versal parameter, that regulates the transition and is estimated 
by means of a number of comparisons with observations, with a 
value ao s 1.2x 10~ 10 m/s 2 . However, the theory does not say 
anything about the form of the function fi(x), which is highly 
degenerate. 
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In our case, we have that the asymptotically constant periph- 
eral velocity v,x, can be inserted into the formula for the cen- 
tripetal acceleration 



P 

so that using our expression in Eq. |20]for v^, we obtain 

2GM 



lap - Ac 



M 



\ 1O 1() M 



P 



(45) 



(46) 



from which we can extract M and insert it into the expression 
of p = p(a). (Of course, the MOND force should now be inter- 
preted as an effective one.) 

There are two solutions, and we therefore obtain the func- 
tions: 



Here, m is the mass of the star, r s is a scale radius, and p(m) is 
a function of the mass to be determined. Since we consider all 
stars as being equal, p is "universal" and depends on the mass 
only. 

The situation is very different for a spiral galaxy, first, be- 
cause it is not spherical and second because it is not a continu- 
ous distribution of mass. Thus, the correct procedure would be to 
start from a metric with cylindrical symmetry, and to solve the 
resulting fourth order equations. A very rough attempt would 
consist of computing the force acting on a test particle near the 
edge of the disk, by considering the galaxy as made up of, say, 
10 10 stars like the Sun, (without either bulge or intergalactic 
dust, which is indeed very crude), and adding the forces. 

The force acting on a unit mass test particle is therefore 



2„4 



A l c 



a = 



4G x 10 10 M o ' 
From the MOND estimate of ao we then find 
A = 2.8- 1(T 7 , 



(48) 



(49) 



in very good agreement with the result extracted from the above 
fit. 

The great advantage of this formulation is that we obtain the 
universal function: 



p(x) = 
which can be tested. 



1 + 2x - Vl +4x 

2x 



(50) 



6. Point-like and extended bodies 

As discussed in previous sections, considering the corrections 
to the potential from galaxies to Solar System scales, with the 
same values of parameters, is rather unsatisfactory. Here we 
want to discuss possible explanations and take into account also 
extended bodies. 

First, both the Sun and a galaxy are not point-like sources. 
In the case of the Sun, however, assuming spherical symmetry 
and neglecting as usual the rotational contribution, we may still 
use the rotational invariant form of the metric in Eq. ©. There 
are important differences if the source is not point-like: other el- 
ements may enter in the computation of the integration constant, 
for instance the radius of the object. The correct treatment would 
be to solve for the internal metric and join it smoothly with the 
external one, so as to take care of the change in density with 
radius. A reasonable approximation would be to assume a uni- 
form density, so that only radius could be involved. Moreover, 
we may consider for the moment all stars as being equal to the 
Sun, so that the correction to the potential can be expressed again 
in terms of the mass 



Gm 

<£ = + p(m) log(p/r s ) . 

P 



(51) 



F = 



Gm p(m) 



(52) 



8aG x 1O 1() M + A 2 c 4 ± Ac 2 ^\6aG x 1O 1O M + A 2 c 4 

u±(a/an) — 77 x 

^ ' ' 8aG x 10 10 M Q 

(47) 

A first important remark is that p is indeed independent of the 
source mass as well as of distance, and it is therefore universal, 
as requested. By comparison with the asymptotic limits, on the 
other hand, we have to choose p_ and get 



Let the test particle be situated on the x-axis, for instance 
(and of course on the galactic plane), R be the radius of the 
disk, and N the number of stars. The first term gives the usual 
Newtonian force, so that we only have to sum the correction. 
Due to the non linearity of the dependence on the mass, this is 
not a straightforward task. We give here just a rough estimate 
(a more precise computation can be made and gives the same 
answer). 

Because of the symmetry, the total force is clearly directed 
towards the center, so that we need to compute only the x com- 
ponent. The correction is 



Zp(m) 
l —— cos(0;) , 
, Pi 



(53) 



but, since < p, < R and < cos(#,) < 1, we may substitute a 
mean value and write 



8F 



T P-(m) 
R 



(54) 



where ^ is a number of order 1 . From the above discussion, on 
the other hand, we have seen that a good empirical expression 
for a galaxy is given by 



5F : 



R 



(55) 



where A* is the same as the A used before, but in appropriate 
units. Comparing the two expressions we get 



p(m) 



X N 



A, yfm 



(56) 



X 



Remembering that x ~ 1 , we see that the force law for the 
Sun is the same that we used for a galaxy, but with a coefficient 
lowered by a factor ~ Vn. In other words, the suppression of 
the correction acts in a twofold manner, because of the very low 
mass of the Sun compared to a galaxy, and because it is a sin- 
gle "almost point-like" object, instead of a compound one. The 
consequence is that all the numerical values of the local tests 
presented above should be lowered by a factor ~ 1CT 5 . In this 
case, the approach works very well. 

The main problem of this argument lies not so much in the 
rough estimate of x, but in the fact that the number of stars is not 
the same for all galaxies, and that the stars may be very different 
in mass and size. A more satisfactory computation should be 
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done starting from Eq.d52l), with a reasonable guess for fJ.(m) 
(may be fi(m, r, ...)), and computing the rotation curves with the 
right distribution of luminous matter and number of stars. This 
is not a simple task and is postponed to future work. 

7. Concluding remarks 

Starting from a reasonably simple f(R) model for gravity, we 
have shown that it is possible to obtain promising astrophysical 
results, which do not require dark matter. What is most impor- 
tant is the fact that, with the same model, it appears possible to 
do the same also at a cosmological level, where exact s olution s 
have been discussed preliminarily in Capozziello et al. (|2008b). 
Of course, much work has still to be done on both scales. The 
work here is only indicative of the concrete possibility of test- 
ing a specific f(R) model of gravity on astrophysical grounds. 
For cosmology, at least the structure formation and the cosmo- 
logical microwave background radiation spectrum must still be 
investigated. On local astrophysical scales, on the other hand, 
more realistic models for objects like galaxies, for example, are 
necessary. 
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